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Abstract 

The 2D gravity described by the action which is an arbitrary function of the 
scalar curvature f(R) is considered. The classical vacuum solutions are analyzed. 
The one-loop renormalizability is studied. For the function / = RlnR the model 
coupled with scalar (conformal) matter is exactly integrated and is shown to de- 
scribe a black hole of the 2D dilaton gravity type. The influence of one-loop quan- 
tum corrections on a classical black hole configuration is studied by including the 
Liouville-Polyakov term. The resulting model turns out to be exactly solvable. The 
general solution is analyzed and shown to be free from the space-time singularities 
for a certain number of scalar fields. 
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1 Introduction 



Recently, much attention has been paid to the investigation of models of two- 
dimensional (2D) gravity. It is well known that the Einstein-Hilbert action in two 
dimensions coincides with the topological Euler number and, therefore, does not 
determine any dynamics for gravitational (metrical) degrees of freedom. Hence, one 
should consider some alternative dynamical descriptions of 2D gravity. One of the 
simplest models, mainly inspired by the string theory, is dilaton gravity [1], gravita- 
tional variables are the dilaton and metric fields (4>, g^ u )- In empty (without matter) 
space the classical equations of motion are exactly integrated [1] and the solution 
describes the 2D black hole. On the quantum level, it has been shown that this 
model is renormalizable [2]. The coupling with conformal matter is again exactly 
solvable classically and the solutions are configurations describing the formation of 
a black hole by collapsing matter [3] . 

An other way is to formulate the theory of 2D gravity in the framework of a 
consistent gauge approach. Independent variables are now vielbeins and the Lorentz 
connection (e a ,uj a b ). The theory with Lagrangian quadratic in curvature R and 
torsion T [4] was shown to be exactly solvable [5]. One class of the solutions contains 
the de Sitter space-time with zero torsion. Other solutions are of the black hole type 
[5]. Generally, one can consider the Lagrangian to be an arbitrary (not-necessarily 
quadratic) function of curvature and torsion [6] . Such a theory has essentially the 
same type of classical solutions. 

Describing the gravitational degrees of freedom on the 2D manifold M 2 only 
by the metric (g^ u ) without introducing any additional variables, one considers the 
following action: 



where f(R) is, in principle, arbitrary (non- linear) function of the scalar curvature R 
determined with respect to the 2D metric g^ v . Theories of such type were studied 
in higher dimensions [7] and in two dimensions [8,9]. Was observed that the theory 
(1.1) is equivalent to some type of scalar-tensor ((f), g^u) theory of gravity. Moreover, 
it was shown in [10] that (1.1) with Lagrangian f = RlnR describes the same black 
hole space-time as the string inspired 2D dilaton gravity. 

One of the motivations for recent investigations of 2D gravity (mainly of the 
dilaton type) that it can be considered as a "toy" model to study the process 
of formation and subsequent evaporation of a black hole. It has been argued by 
Hawking [11] that such a process is not governed by the usual laws of quantum 
mechanics: rather, pure states evolve into mixed states. However, it is commonly 
believed that a successful quantization of gravity and matter will provide us with a 
consistent solution of this problem. Quantum corrections may completely change the 
gravitational equations and the corresponding space-time geometry at the Planck 
scales. This problem is hard to analyze in four space-time dimensions. However, in 
two dimensions one can attempt to attack this problem using the dilaton gravity 




(1.1) 
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theories as a toy model [3] . These toy models have an explicit semiclassical treatment 
of the back reaction of the Hawking radiation on the geometry of an evaporating 
black hole by including the one-loop Polyakov-Liouville term in the action (the 
review can be found in [12]). Unfortunately, the resulting equations are not exactly 
integrated and one can not obtain a definite answer. Therefore, one can try to 
find another theory of 2D gravity (among the alternatives) for which the relevant 
semiclassical equations would be analytically solvable. 

The main goal of our paper is the study of this problem for 2D gravity described 
by an action of the form (1.1) along the lines of ref.[3]. We show that for f = RlnR 
the semiclassical field equations are exactly integrated and one can obtain a definite 
answer about the structure of space-time when the backreaction of the Hawking 
radiation on the black hole geometry is taken into account. 

This paper is organized as follows. In the next two sections we investigate some 
aspects common for theories described by the action (1.1). In Sec. 2 we demonstrate 
the integrability of classical field equations and find the exact solution. The one- 
loop renormalizability of the theory is analyzed in Sec. 3. In the next two sections we 
mainly consider the case / = RlnR. The coupling with conformal (scalar) matter 
is shown to be exactly solvable classically in Sec. 4. The backreaction is taken into 
account in Sec. 5. 

2 Classical solution of the model 

Under variation of the action (1.1) with respect to the metric g^ v we obtain the 
following equations of motion: 

v M v,[/'] = \g^{f(R) - Rf'(R) + 2D(/')} (2.1) 

where /' = d R f(R) and □ = V%. 

At first sight, (2.1) is a system of differential equations of very high order with 
respect to derivatives. For example, if / = R 2 , then (2.1) are equations of fourth 
order of metric g^ derivatives. However, we will see that it is not really so and the 
system (2.1) is rather easily solved. 

Let us analyze at first possible solutions of (2.1) with the constant curvature 
R = R\ = const. In this case we obtain that f'(R) = const = f'\n=R 1 everywhere 
in M 2 . Then, from (2.1) we get that such a solution exists if the function 

V(R) = f(R) - Rf'(R) (2.2) 

is zero at the point R = Ry. V(R\) = 0. If V(R) becomes zero at P different 
points Ri,i = 1,2, ...,P, then for given f(R) there are P different solutions of (2.1) 
with constant curvature. An additional condition is that the function f'(R) must 
be finite at R = R{. 

Assuming R to be a non-constant function on M 2 , we consider a new variable 
4> = f'(R) provided that this equation is solved (at least locally) with respect to R: 
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R = R{4>). Denote V(<p) = V(R((f>)). Then (2.1) is rewritten as an equation on the 
new field 4>: 

V„V„0 = \g^{V(4>) + (2.3) 

We obtain from (2.3) that £ M = e^d u <p is the Killing vector [13]. Consequently, 
the field cf) can be chosen as one of the coordinates on M 2 . Then, metric reads 

ds 2 = g(cf>)dt 2 - -±^dxj> 2 (2.4) 

From (2.3) we get that 

n<f> = -V(<f>) (2.5) 
For the metric (2.4) we have □</> = —g'(4>) and eq.(2.5) reads 

d^tt) = V{tj>) (2.6) 

The solution takes the form 

g(4>) = A + J VtyW (2.7) 

One can see that our model (1.1), (2.1) seems to be equivalent to some kind of 
2D dilaton gravity with the dilaton field (f> and potential V((f>). 

Thus, surprisingly, our initial higher-derivative equations reduced to the first 
order equation (2.6) independently of the concrete form of the function f(R). As a 
result, the solution (2.7) is determined only by one arbitrary integrating constant A. 
The Killing vector ^ = e^d v (f> has bifurcation at a point where £ 2 = — (V</>) 2 equals 
zero. One can see from (2.4) that g(<f>) = at this point and we have a horizon. 

Since the scalar curvature for the metric (2.4) is equal to R = —</'(</>), one can 
easily check that curvature for the solution (2.7) really coincides (if f"(R) / 0) with 
R = R{4>) obtained by solving the equation <p = f'(R). 

In the vicinity of points where <fi' R = f"(R) = the equation <j) = f'(R) cannot 
be solved in a unique way. It is the only place where above solution is in-correct. 
Near a point like that the function <j>(R) is as shown in Fig.l : there are two values 
R a , Rb which correspond to the same value 4>. So there are two branches of solution 
of equation <fi = f'(R). Let us consider this in more detail. Let <j>'(R) = for some 
finite R = Ro. Note that only zero of odd order is interesting for us. In the vicinity 
of Ro the function <p'{R) can be represented as follows 

4>'{R) = a(R - Ro) 21 *- 1 , k= 1,2, ... , and a > (2.8) 

and hence we obtain 

HR) = ^(R~Ro) 2k +b (2.9) 

and 

^ = miTT ) (R - lt ^ +1 + bR + c <2 - 10) 
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There are two branches of the solution of eq.(2.9) with respect to R: 

R = R o ±[(0-b)—]2k (2.11) 
a 

In the vicinity of the point x £ M 2 where R(x) = Ro there are two regions: where 
R > Ro and where R < Ro- Our solution (2.4), (2.7) is valid in any of these regions 
taken separately. Consider the region where R > Ro. Suppose for simplicity that 
6 = 0, a > 0; then <p > 0. Then, we get for the potential V: 

V&) = -^ T (-)^0 1+1/2fe -R^ + c (2.12) 
2k + 1 a 

The corresponding metric (2.7) for cf) > reads 

m = a - + * - g^ ^r* (2.i3) 

One can see that the metrical function #(</>) (2.13) has regular in <f> = first and 
second derivatives 



</(0) = 0, /(0) = -#0 

However, the non-analyticity of (2.13) in = manifests itself in that all the 
following derivatives are singular at this point: 

g {p) (0) = ±oo, p > 2 

This singularity means, in particular, that invariant (Vi?) 2 is singular at the 
point x where R(x) = Ro- It should be noted that singularities of this type were 
earlier observed in [14,15] for rather different theories. 

In the region where R < Ro we get 

R = Ro - {^f/ 2k 

v{(t>) = 2TTI ( T )1/2V+1/2fe " Ro(t> + c ' ^ > (2 ' 14) 

Thus, the total space-time in the vicinity of the point R = Rq (<j> = 0) is repre- 
sented by gluing of two sheets (the coordinate eft > can be used to parameterize 
the points of both sheets in the neighborhood of <j> = 0). The total space-time is 
shown in Fig. 2. 

Really, the scalar curvature R itself can be used as one of the coordinates. It 
covers, in particular, the whole vicinity of the point Rq- Then, in the coordinates 
(t, R) the metric reads: 

ds 2 = g(R)dt 2 - lf "™ 2 dR 2 (2.15) 
9{R) 
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For R ~ Rq we can put g(r) ~ 1, f"(R) ~ a(R — Ro) 2k and hence the metric takes 
the form 

ds 2 = dt 2 -a 2 {R- R ) 4k dR 2 (2.16) 

Let us now consider some examples. 
Example 1. 

f(R) = R]nR (2.17) 

In this case <j) = hiR + 1, R = e^ -1 . Hence 1/(0) = -R = -e^ -1 . Since the 
potential V(-R) = — i? is zero in i? = 0, it seems that one of the solutions is flat 
space-time. However, the function /'(i?) = lnR+ 1 is not defined for i? = 0. So if 
we come back to eq.(2.1), we observe that flat space-time is not really a solution of 
the field equations. 

If R is a non-constant function on M 2 , the solution is given by the metric (2.4) 
with 

g{4>) = A - R(4>) = A - e* -1 (2.18) 

This solution coincides with that obtained in 2d dilaton gravity and describes asymp- 
totically flat black hole space-time. The essential difference of the solution (2.18) 
from that we have in dilaton gravity is that it doesn't describe flat space-time for 
any integrating constant A. 

The Lagrangian (2.17) seems to be ill-defined at R = 0. However, we see that 
curvature R{4>) is everywhere positive and the point R = really lies at the spatial 
infinity. 

Example 2. 

f(R) = aR 2 + bR + c (2.19) 

In this case <fi = 2aR + b, R = ^(4> — b), i.e. R(4>) is linear function. Then, we 
get V(R) = -aR 2 + c = -^{<j> - b) 2 + c. If c/a > 0, then V(R) is zero at the 
points R = ±y/c/a. Thus, there are two solutions with the constant curvature: 
R = ±y / c/a. If R is non-constant on M 2 , then the solution is given by (2.4) with 
g(4>) in following form: 

g ( (j) ) = X + c ( f>--±-( ( f>-b) 3 (2.20) 

This function has extremums at the points <p\^ = b ± 2a^Jc/a corresponding to the 
curvature R = ±y/c/a. Depending on the constant A (if a,b,c are fixed), g(cj)) can 
have one, two or three zeros. It is worth noting that the space-time described by the 
metric (2.20) is not asymptotically flat. R = is reached at the point 4> = b which 
stays on finite distance from any point 4> / ±oo. Thus, the points (f> = ±oo lies at 
asymptotical infinity and the curvature is singular at this point. In this sense, the 
solution (2.20) is similar to that obtained in the 2D theory of gravity with torsion 
described by the action quadratic in curvature and torsion [5]. 
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It should be noted that a behavior like that is rather typical of polynomial 
gravity (1.1). Indeed, let the function f(R) near R = look like as f(R) = aR a + c. 
Then, = aaR a ~ l , R(<f>) = (^) 1/q_1 . One can see that for a > 1 the point R = 
corresponds to (ft = and, consequently, lies at a finite distance from any point 
0/0. It means that space-time is not asymptotically flat. The last is reached 
only if a < 1: then R — > means — > +oo. However, in this case the function 
/(i?) is not analytical in R = 0. It is the case in Example 1., where the solution is 
asymptotically flat. More generally, the solution is asymptotically flat if the function 
f(R) satisfies the condition: f'(R) — > ±oo if R — > 0. It is easy to see, however, that 
flat space-time is not a solution of field equations in this case, as we have seen in 
Examplel . 



The complete quantization of the model (1.1) is a rather difficult problem. In this 
section, we just calculate one-loop counter-terms and check the renormalizability 
of the model in one loop and not considering these problems as the unitarity. We 
assume in this section that f"(R) ^ 0. 

We use the background method. The metric g^ u is written in the form: g^ v = 
g^iv+hfxv, where g^ is a classical background metric, is a small quantum field. In 
the conformal gauge we have h^ u = a/2g^ v and the theory reduces to quantization 
of only a conformal mode a. Expanding the action (1.1) in powers of a we obtain 
the quadratic in a expression 



where the curvature R = R[g] and the Laplacian □ = -^d^yfgg^ d u ] are deter- 
mined with respect to the background metric g^ v . We see from (3.1) that (1/ /") is 
effectively the loop expansion parameter for gravity. 
The action S q [a] can be written in the form 



3 One-loop renormalization 



S[g^] = S d [g ilv \ + S q [a\ 

SM = J [f"(R)(Oa) 2 - 2Rf"(aaa) + (R 2 f + / - Rf')a 2 ]^d 2 z (3.1) 




(3.2) 



where D is the fourth-order differential operator 



d = (□ + *)/"(□ + y) 



(3.3) 



and the functions X, and Y satisfy the following equations: 



X + Y = -2R, XY = R 2 f" + f - Rf 



(3.4) 
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Calculating the functional integral over the conformal factor a we can compute the 
infinite part of the one-loop effective action 

4c = ^(lncfetD)oo - (lndetA 9/i )oo (3.5) 

where A g h is the standard ghost operator corresponding to the conformal gauge. 

By definition, for an elliptic 2r order differential operator A defined on a two- 
dimensional Riemannian manifold M 2 we get: 

+00 

IndetA = - J jTre~ tA , e +0 (3.6) 

e 

The infinite part is given by (L — > 00): 



r 1? 

(IndetA)^ = -(B L 2 + 2B 1 L + -B 2 ln — ) (3.7) 



where 



(Tre-* A )^ = E ^* (fc " 2)/2r + 0(Vt), 

k=0 

B k = J b k (A)^d 2 z+ J c k (A)^dr, (3.8) 



M 2 dM 2 

b k (A) are the Seeley coefficients for the operator A (&2p+i = 0). For simplicity we 
will assume that M 2 is a manifold without a boundary (dM 2 = 0) and will neglect 
all boundary effects. 

Note that L 2 and L dependent terms are automatically absent in the dimensional 
or (-function regularization. So only the last term in (3.7) is of interest for us. 

Now consider the Seeley coefficient for the elliptic fourth-order operator D (3.3) 
(see [16]). Suppose that for some A4 

A 4 = A 2 A' 2 , detA 4 = detA 2 detA' 2 (3.9) 

Then, we get the corresponding expression for the Seeley coefficients [16]: 

2B 2 (A 4 ) = B 2 (A 2 ) + B 2 (A' 2 ), (3.10) 

Since the operator D (3.3) has the structure (3.9) we obtain that 

2B 2 (D) = B 2 (a + X) + B 2 (f"(a + Y)) (3.11) 

For the operator 

A 2 = a + X (3.12) 
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the following result is well-known: 

6 (A 2 ) = -L ; 6 2 (A 2 ) = ^-{1/6R + X) (3.13) 

To calculate the Seeley coefficients for the operator 

A' 2 = /"(n + F), (3.14) 

it is useful to observe that this operator can be transformed to (3.12) by introducing 
a new metric ~g PhV = (/") _1 £^: 

A' 2 = a~ g + f"Y, (3.15) 
where = -^d^-^/gg^ d v \. Now using the result (3.13) we obtain 

B 2 (A' 2 ) = -L y (1/6R + f"Y)^gd 2 z (3.16) 

M 2 

or in terms of the old metric g^ u eq.(3.16) takes the form 

B 2 (A' 2 ) = i- J (1/6R + Y)^d 2 z (3.17) 

M 2 

Thus, we obtain for B 2 (D): 

^2p) = ^ / (1/3-R + X + Y)^d 2 z (3.18) 

M 2 

Using (3.4) we finally get 

B 2 (D) = ~ f 5/6R^d 2 z (3.19) 

M 2 

The Seeley coefficient for the ghost operator A g h is well-known [17] 

& 2(A^) = ^(2/3i?) (3.20) 

Taking into account (3.7), (3.19-20) we obtain from (3.5) that the corresponding 
one-loop counter-term 



a J R^g~d 2 z (3.21) 



let 

M 2 

is surprisingly non-dependent on the concrete function f(R). So the model (1.1) 
seems to be one-loop renormalizable. Of course, if one takes the next loops this 
result could be changed and possibly not for any f(R) the theory is renormalizable. 
However, we do not consider higher loops here. 
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4 Coupling with conformal matter 

Let us consider interaction of higher derivative gravity (1.1) with 2d conformal 
matter described by the action 

Smat = J \{Vi>?Vgd 2 Z (4.1) 
M 2 

Then, we get the complete system of equations of motion: 
Tfu, = V„V„0 - \g„AV{<i>) + 2D^] + \{d^d v ^ - ^Md a 4>) = (4.2) 
where <f> = f'(R). The equation of motion for matter reads 

□V> = (4.3) 

We will use the conformal gauge in which the components of metric: g ++ = g = 

0; g + _ = \e° '. Then, eq.(4.2) takes the form 

T±± = d±d±<p - d±ad±(t> + -d±Y><9±Y> = (4.4) 

T+_ = <=> 4d+d-(f) = -V{(t))e a (4.5) 
Moreover, we have the self-consistency condition : 

Ad+d-a = R{cp)e a (4.6) 

Equation (4.3) takes the form 
and the solution reads 

ip = ip+{x + ) + i/j^(x~) (4.7) 

For the function f(R) of general form eqs. (4.4-6) are extremely non-linear differ- 
ential equations which are not exactly solved in general. However, in some particular 
cases, for concrete f(R), this problem can be essentially simplified. 

Let us consider the case when the following equation is valid: 

d + d-{a - (f>) = (4.8) 
It is the case when f(R) satisfies the equation: 

R + V{R) = R + f - Rf = 0, (4.9) 
i.e., when f(R) = RlnR. This case was considered in Example 1 . 
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Then, we get 

a - <j) = w + (x + ) + w-(x-) (4.10) 
On the other hand, one can see that (a + (f>) satisfies the Liouville equation 

d+d-{a + <l>) = ±f+* (4.11) 
which has following general solution: 

° + * = ln 1 / 1 (X+ 1?P12 ) = fl* + > *")> ( 4 - 12 ) 

where A and £> are still unknown functions of x + and x~~ respectively. 
Thus, we obtain: 

a = ±p+l w , 0=1/3-^ (4.13) 

From this we see that 

d 2 + 4> - d + ad + <f> = \{d\H - l/2(d + (3) 2 ) - \{d 2 + w - l/2(d + w) 2 ) 
On the other hand, one can see the following identity: 

(d 2 + P-l/2(d + P) 2 ) = {A-x + } 
where we introduced the Schwarzian derivative 

Then, we get for the (++)-component of equation (4.4): 

{A - x + } _ {d l w+ _ l/2(d + w + ) 2 ) + 2Tt + = (4.15) 
where I+ + = l/2(d+V>+) 2 - Similarly, we obtain for the ( )-component of eq.(4.4): 

{B; x~} - {d 2 w- - l/2(d_w„) 2 ) + 2T^_ = (4.16) 

Using the known property of the Schwarzian derivative (see for example [18]), one 
can see that eqs.(4.15)-(4.16) are invariant under SL(2, R) © SL(2, R) group trans- 
formations: 

aA + b 

A — > — , ad — be = 1 

cA + d 

^ mB + n 

B — > , mp — kn = 1 (4-17) 

kB +p 

Under the coordinate transformations x ± — > 7/ ± (x ± ) we have: 

/3(x+,x-) - /?(y + ,y-) - (4-18) 
11 



On the other hand, w± transforms as follows: 

w ± (x ± )^w ± (y ± )-\n(^) (4.19) 

We can use this symmetry to put w± = 0. Then, one obtains equations on functions 
A and B: 

{A-x + } = -{d^ + f 

{B;x~} = -(d-i)^) 2 (4.20) 
When the matter is absent (T±± = 0), the solution of equations 

{A;x + } = 0, {B;x-} =0 (4.21) 
is one of the following types. If A" = 0, then 

A = ax + + b (4.22) 

if A" ^ 0, then 

A = f - 

x + + b 

Correspondingly, we get for B(x~): 

B = mx~ + n (4.24) 

or 

B = d- 

x + n 



A = f-^7T ( 4 - 23 ) 



B = d- m (4.25) 



The metric takes the form: 



r a/ n'U/2 

ds 2 = \ jr^-dx+dx- (4.26) 



(1 



AB 
4e 



Shifting x + ,x~ on constants we get b = n = in (4.22-25). Though A, B depend 
on the set of constants, the metric (4.26) depends only on one arbitrary constant. 
Let, for example, A" = B" = 0, then 

, 9 cdx + dx~ , , . 

ds = 7] c2 + ; ( 4 - 27 ) 

where c = ^am. In other cases, if A" and B" are not zero, the metric takes the 
form: 

<fe 2 = -— r (4.28) 

where c = (1 — ^)(am)~'/ 2 , The scalar curvature is given by the formula: 

R ( 4 29) 

R ~ e (1 - 4f ) ( 9) 
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It has singularity if AB = 4e, which one can also see from eqs.(4.27), (4.28). The 
points of horizon satisfy: AB = 0. The space-time (4.27-28) is of the same type as 
the black hole solution in the 2D dilaton gravity [1,3]. However, there is no such 
integrating constant for which the metric (4.27-28) is flat. The flat space-time is not 
a solution of field equations that has already been noted above. This is an essential 
difference between the string inspired 2D dilaton gravity [1] and higher derivative 
gravity (1.1) with / = RlnR. Therefore, eqs.(4.4)-(4.6) do not describe the black 
hole formation from regular (flat) space-time due to the infalling matter as we had 
in the dilaton gravity [3]. The "bare" black hole is necessary. The infalling matter 
only deforms this initially singular space-time. 

As an example let us now consider the falling of 5-like impulse of matter on the 
black hole. The matter energy-momentum tensor takes the form: Tt, = XS(x + — 
Xq), = (A > 0). It describes the 5-like impulse of matter propagating along 
the ^"-direction. Suppose that the space-time for x + < Xq is a solution of the field 
equations without matter such that A = ax + , B = mx~ . For x + > Xq the function 
B{x~) is the same while A(x + ) is found from the equation: 

{A; x + } = -X5(x + - x+) (4.30) 

For x + > Xq the function A(x + ) is a solution of eq.(4.30) with the zero right-hand 
side 

A = aX ' +l3 , arf-K0 = l (4.31) 
KX + 7 

where the constants a, (3, k, 7 are found from the continuity condition of functions 
A(x + ), A'(x + ) and the gap condition for A"(x + ) at the point x + = Xq. The last 
condition is easily obtained integrating (4.30) in the interval (xq — e, x^ + e) and 
then taking the limit e —> 0. As a result one obtains: 

A"{xt + 0) - A"(x^ - 0) = -\A'(x+) (4.32) 

From continuity of A(x + ) and A'(x + ) one gets: 

aa?o = — ^— (4.33) 
KXq + 7 

a = ( KX + + 7 )-2 (4.34) 
and from the gap condition (4.32) we obtain 

J^TT)- X (435) 

These equations and a~f — K(3 = 1 are enough to find the form of A{x + ) for x + > Xq: 



A(x+) = a X+ + f { f~^ (4.36) 
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The metric for x + < Xq takes the form 



, o Jamdx + dx , , . 

* 2 = ( i_ ^ x+I - } 



and the corresponding curvature is the following 



„ i J am 

*-«-' (1 -W> (438) 



This metric describes the black hole with horizon in x + x =0 and singularity at 

4- — 4p 
2V X = — . 

am 

We will assume that Xq > 0, i.e. impulse falls from asymptotically flat region 
which lies right of horizon (x + = 0). Then, for x + > Xq we obtain for the metric 

, 9 Jam dx + dx~ 

ds 2 = ^ t 4.39 

and the curvature 

= e^y^l + ^(x+ - ar+) - ^(x+ + ^(x+ - x+))x-]- L (4.40) 
One can see from (4.40) that for x + > Xq singularity lies on the curve: 

x - 4e (i + 4)) (441) 
«™ (a;+ + ^i( x +_ x + ) ) 

The derivative of the function (4.41): 

5 + x- = -— (x+ + ^(x+-x+)r 2 

is negative and we have that for x + > Xq the function (4.41) is the monotonically 
decreasing one smoothly glued with x~ = ^ at x + = Xq. Moreover, in the limit 
x + — > oo it limits to x~ — > x^ = ^(xq~ + 2/A)^ 1 . The total space-time for all x + 
is shown in Fig.3. In the asymptotically flat region (x + > 0) it is similar to that we 
have for the 2D dilaton gravity case [3]. 

It should be noted that the function f(R) = RlnR is not a unique one for which 
equations (4.4-6) are exactly integrated. Indeed, we obtain from (4.5-6): 

4d + d-(a-(p) = {R + V)e a 

4d+d- (a + <f>) = (R- V)e a (4.42) 
These equations are reduced to the system of the Liouville equations if 

R + V = cuT*, R-V = be^ (4.43) 
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where a, b are constants. These conditions are equivalent to the system of differential 
equations on the function f(R) 

R + f-Rf' = ae~ f ' 

R- f + Rf = be f ' (4.44) 
One obtains immediately from this 

R= l/2[ae- f ' +be f '} (4.45) 
Solving (4.45) with respect to /' one obtains: 



/- = ln g±V f^ (4.46) 



Integrating this we finally get: 



R ± VR 2 - ab 



f{R) = R\n ^ — T VR 2 - ab (4.47) 

where both signs (±) are available if (ab) > 0. We do not consider here this kind of 
theory. Note only that for a, b ^ it describes the asymptotically singular rather 
than asymptotically flat space-time. 



5 Solution with backreaction 

As we have discribed in the Introduction quantum corrections are usually assumed 
to remove the black hole singularity One can try to analyze this problem semi- 
classically considering quantum gravity coupled to a large number N of free scalar 
fields. In the limit fi — > with NU fixed, it is a system in which the leading order 
of a perturbative expansion is a quantum theory of matter in classical geometry. 
Integrating out the matter we have an effective action describing the backreaction 
of matter and Hawking radiation on the geometry, which we hope to treat classi- 
cally [19]. In ref.[3] it was proposed to use this approach to study the problem in 
two dimensions for dilaton gravity. However, the resulting quantum-corrected field 
equations are not exactly solved [3,12,20] though some reasons observed in favor 
of that singularity are still present in this semiclassical theory. We apply here the 
approach of [3] for theory of gravity described by the action (1.1). 

In two dimensions, integrating out the conformal scalar fields one gets the 
Polyakov-Liouville action: 

S pl = iL i d 2 x ly f^ ( d 2 x 2 V=gR(x 1 )n- 1 (x 1 ,X2)R(x 2 ) (5.1) 
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here □ denotes the Green function for the Laplacian. It should be noted that 
Spl incorporates both the Hawking radiation and the effects of its backreaction on 
the geometry. We neglect here the contribution of the ghosts [21]. The full effective 
action 

^eff = S gr + S ma t + Spl (5.2) 

gives rise to the following system of equations (the metric is taken to be conformally 
flat = \e a ): 

d ± d±(f> - d±ad±<f> + 2c[d 2 ± a - ^(d±af - t±] + T| ± = (5.3) 



Ad + d.(k = -e a (V{(t)) + 2cR) (5.4) 

where c = || and T±± = ^d±tpd±tp. Equation (5.4) is obtained as variation of the 
action (5.2) with respect to </+_. Since the scalar curvature is a known function of 
4>, we must add the condition of self-consistency: 



The scalar matter equation 



gives 



4d + d-a = R{<j))e a (5.5) 



For general function f(R) these equations seem to be not exactly integrated. 
Therefore, we will consider in this section only the case 



S gr = J R\nR^f^)d 2 x 



and show that for this type of gravitational action the system (5.3-5) is exactly 
solved. In this case R(<f>) = e^ -1 , V(<j>) = —R. Equations (5.4), (5.5) take the form 



(5.6) 



Ad + d-a = (5.7) 



Let c / 1. 

Then, from (5.6-7) we obtain 



d + d-[(l - 2c)a - (/)] =0 

d+d„ [</> + a] = iiZ-^e^ ( 5 . 
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These equations are easily solved 



(1 — 2c)a — <j) = w + (x + ) + w- (x ) = w 
A'B' 



[<f> + a]=ln - — = P ( 5 - 9 ) 

4e 



(1 - ^ABf 



where A = A(x + ), B = B(x ). Finally, we get for the conformal factor a and field 
0: 

* = %hr^ + = tr^y" - %hrf < 5 - 10) 

With respect to the coordinate changing x ± — > y ± (x ± ) we have 

/3(x + ,x") -> ln/?(y + ,jT) + (d+y + d-y~) 

w ± (x ± ) -^w ± (y ± ) - (l-2c)ln(d ± y ± ) (5.11) 

Hence, the fields <r and transform as a usual conformal factor and a scalar field, 
respectively: 

a(x + ,x~) -> a(y~,y + ) - ln(d^y~d + y + ) , 4>(x + , x~ ) ^ (f>(y + , y~ ) (5.12) 
One can easily see that 

S|0 - d±ad ±< /> + 2c[4o- - ^(9±a) 2 ] = 

1 [dlP-hd ± P) 2 ] + —^—-[d 2 ± w(-l + 2c) + kd ± wf] (5.13) 



2(1 -c) L * 2 V ^ 7 J 2(1 -c) L * v ' 2 

As before, we have in terms of the Schwarzian derivative 

dlP-\(d + pf = {A-x + } 

d 2 _f3-\{d-(3f = {B-x-} (5.14) 

Let moreover c / 1/2, then we can use the symmetry (5.11) to put w = 0. Then, 
equations (5.3) take the form 

{A;x+} = -2(1 - c)Tt + + 4c(l - c)t + (x+) 

{B;x-} = -2(1 -c)Tt_ +4c(l - c)*_(0 (5.15) 

Eqs.(5.15) are ordinary differential equations with respect to A, B. 
The metric and curvature, respectively, read: 

A 1 B f 1 
ds 2 = e^^dx+dx- = \ W^dx + dx- (5.16) 
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1 (l-2c)/3 1 A>B> (l-2e) 

R = -e^^V = -[ — l^y (5.i7) 

If the right-hand side of eq.(5.15) is zero (i.e., matter is absent), then the solution 
of eq.(5.15) is already known (see (4.22-25)). Let, for example, it take the form 
(4.22), (4.24) with b = n = 0: A = ax + , B = mx~ . Then, the metric and scalar 
curvature take the form 

ds> = [—j^—^dx^dx- (5.18) 

R=-[ TT^ (5.19) 

If c = 0, we obtain the "classical" black hole space-time with space- like singularity 
at x + x~ = — (we assume that am > 0). 

am v / 

For c < 1/2 or c > 1 we see from (5.19) that space-time still has singularity at 
x + x~ = am ^i_ c ^ ■ It is time-like for c > 1 and space-like for c < 1/2. As before, 
the points of horizon satisfy the condition (dcf)) 2 = 0, which for (5.18), (5.19) means 
that x + x~ = 0. The diagram of this space-time is shown in Fig.4. The regions I 
and III are asymptotically flat. 

Let c > 1 and consider the falling in this space-time of the matter impulse 
T++ = XS(x + — Xq), A > 0. We will neglect contribution of t±(x ± ) into (5.15). We 
assume that impulse falls in the region I which is asymptotically flat, i.e. Xq < 0. 
For x + < Xq the space-time is described by the metric (5.18) and has curvature 
(5.19). For x + > x^ the solution of eq.(5.15) is found in the same way as before 
(see the previous section). Moreover, the solution has the form similar to (4.36) 

A(x+) = a X + ^ ( ; x- > *+ (5.20) 

We obtain correspondingly for the metric 

ds 2 = H^jl + ^(1 - c)(x + - x%) - 

am( \~ 6) x-{x+ + ^(1 - c)x+(x + - x+))\^dx + dx- (5.21) 
and scalar curvature 

R = l/e(om)l^[l + ^(1 - c){x + - atf) - 

am{ \- C) x-(x+ + \{l - c)z+(z+ - 4W 2 -^ (5-22) 
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For x + > Xq the singularity lies on the curve 

x - = 4e 

am(l-c) (x+ + ^)£i (x+ _ x + )) 

which is smoothly glued with the curve x~ = am fi_ c ) l/ x+ a t the point x + = Xq. 
Calculating a derivative of the function (5.23), we obtain that 

d + x ~ = - „ m n-J x+ + ^T^ x o( x+ - 4)r 2 (5-24) 

i.e., the function (5.23) is monotonically increasing (remember that we consider the 
case c > 1). The function (5.23) takes an infinite value at = (1 + j^i~^^t) x o i 

Xq < x^ < 0. It means that singularity in the region II is slightly shifted for 
x + > Xq , as is shown in Fig.5. On the other hand, the singularity in the region 
IV asymptotically tends to x~ = f^(l + ^"^ Xq The resulting space-time is 
shown in Fig.5. We see that for large N , c > 1, the singularity doesn't disappear 
but simply becomes time-like, which is similar to that we have for the 2D dilaton 
gravity [20]. 

Another case happens if c lies in the interval 1/2 < c < 1. One can see that 
power in the expression for the curvature (5.19) becomes negative. Hence, the 
metric (5.18) describes space-time which is regular for any finite x + and x~ . In 
particular, it is the case for the points on the line x + x~ = am ^i_ c ^ (or AB = jfz^)- 

The curvature is zero though the metric takes an infinite value on this line. 

We obtain singularity if x + or x" takes infinite value. It is convenient to change 
variables: x ± = (y ± ) _1 . Then, the metric and curvature take the form 

7 2 1 r i^- , + , - 

ds = -, — ; rrr 77 — t \ 1 - c dy ay 

« = l/e[^(l-^^)]^ (5.25) 
yam 4ey+y~ 

In the coordinates (y + ,y~) the singularity lies on the light cone y + y~ = 0. Asymp- 
totically (for y + y~ — > oo) this space-time is of constant curvature. 

The special case is c = 1/2. 

One can see from (5.11) that (a; 1 * 1 ) transform as usual scalar fields. Hence, 
one cannot put w = 0. Taking into account (5.13), (5.14) we obtain for eq.(5.3): 

{A-x + } + \(d + wf+Ti + -t + = Q 

{B-x~} + -{d-wf + T^_-i_ = (5.26) 
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The metric and curvature take the form: 

A' B'e w+ e w ~ 

ds2 = n _ ab\2 dx+dx ~ 

y 1 8e I 

R = l/ ee - w+ e- w ~ (5.27) 

We may use the coordinate freedom to choose A and B as new coordinates: u = 
A(x + ), v = B{x~). Since under x + — > y + (x + ) the Schwarzian derivative transforms 
as follows [18]: 

{A-x + }^(d + y + f{A-y + } + {y + -x + } (5.28) 
eqs.(5.26) are rewritten in the following form: 

l -( du w) 2 +Tt = t u 

±(d v w? + Tt = t v (5.29) 

Note that inhomogeneous piece of law (5.28) cancels in (5.26) with the corresponding 
transformation of t±, so that finally we come to expression (5.29). In the new 
coordinates we have 

e w+{u) e w-(v) 

ds2 = n _ m\2 dudv 

R=l/ee- w+ Me- w ~M (5.30) 

In general, the solution of eqs.(5.29) depends on the choice of boundary condi- 
tions, i.e., on appropriate functions t u , t v . These functions mean the flow of the 
Hawking radiation due to the falling matter with energy-momentum tensor Tff v . 
Physically, it seems to be reasonable to assume that the energy back radiated can- 
not be larger than the energy of the falling matter, i.e. t± < r t±±- From (5.29) 
we obtain that unique possibility is the following: T±± = t±. Hence, one gets 
w± = const and, consequently, the total space-time is of the constant curvature: 
R = e~( w+1 \ Notice, that only for c = 1/2 there exists a constant curvature (de 
Sitter) solution of equations (5.3)-(5.5). 

The other special case is c = 1. 

Then, as one can see from (5.6) and (5.7) we obtain 

(a + (/)) = (5.31) 

This equation has the solution 

cr + (f> = w = w + (x + ) +w~(x~) (5.32) 
Inserting this into (5.7) we get an equation on conformal factor a: 

_L e ™ = ^-e w ~e w+ (5.33) 
4e 4e 
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which has the solution: 

v = ^- e j X e w+{z+) dz + J X e w '^dz- + a(x + ) + 0(x~) (5.34) 

We use the coordinate freedom to put a(x + ) = 0, (3(x~) = 0. One can see that 
d±a satisfy the following equation 

d±a = d±w±d±a (5.35) 

Taking this into account and putting (5.32), (5.34) into (5.3) we obtain that w± 
satisfy the following equations: 

d 2 ±w ± = -{Tf± - 2t±) (5.36) 

The general solution of (5.36) takes the form 

du j (T£ ± -2t ± )dz (5.37) 

If matter doesn't contribute (i.e. the right-hand side of (5.36) is zero), then 

w + = ax + + b, w~ = mx~ + d (5.38) 

where a, b, m, d are constants. Below we consider the case a, m > 0. 

Let us now consider the 5-like matter contribution (t± are putted to zero) T^^_ — 
XS(x + — Xq), = (A > 0). Then, w~ = mx~ + d for all x + and w + takes the 
form 

w + = ax + + b, if x + < Xq 

= (a- \)x + + b + \x£, if x + > x% (5.39) 

Choosing the integrating constants in (5.34) to be zero we have correspondingly for 
a: 

a = -J— e ™-+d e ax++b^ if x + < x + 
Aeam 

e mx-+d ][e {a-\)x++b+\x+ _ ^ e ax++b^ x + > x + 



Aem(a — A) 



We see that the metric = ^e a is everywhere positive and regular for any finite 

x + , x~. 

It is worth observing that the scalar curvature R = \/ee w ~ a can be written in 
the form: 

R = -xe~ ax (5.41) 
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where we introduced the function \ (x > 0) taking the form 



X = 




.+ 
o 



(5.42) 



and function a is a 



i 



for x + < xt and 



1 




(5.43) 



a = 



4em(a — A) 



a 



for x + > Xq . One can see that a is positive both for A < a and A > a. Moreover it 
takes positive finite value in the limit A — > a: 



Hence, we obtain that the function a is positive for all x + and the curvature R 
(5.41) is finite for all x + , x~ . We obtain the asymptotically flat space-time which is 
free from singularity and horizons. 

Thus, the solution of equations (5.3)-(5.5) for c = 1 describes everywhere regular 
space-time. This case gives us good example when the quantum corrections (taken 
into account in the form of the Polyakov-Liouville term in the action (5.2)) can really 
remove the space-time singularity of the classical (black hole) solution. It should 
be noted that this result essentially depends on the quantum state or, equivalently, 
on the choice of appropriate boundary conditions (functions t±). In the case under 
consideration the choice was to get asymptotically flat space-time. 

Some remarks are in order. As we have seen in Section 3. (eq.(3.1)) for action 
(1.1) the value (1//") is effectively loop expansion parameter for gravity. The 
semiclassical approximation (5.2) is valid under the condition: |1//"| << N. For 
/ = RlnR we have (f")^ 1 = R. Consequently, we obtain the condition: \R\ << N. 
This condition is rather natural and it means that semiclassical approximation works 
far from the points where the curvature infinitely grows. It is seen from the above 
consideration that for fixed very large N (c » 1) there necessary exists a region 
near the space-time singularity where this condition is not valid and hence the 
semiclassical approximation failed. However, we can see that for N = 48 (c = 1) 
the curvature R (5.41) is bounded and has a maximum value: R m ax = (e 2 a) -1 . We 
have that a = (Aeam)~ l where (am) is an integrating constant. Thus, we obtain 
that a semiclassical approximation is valid for c = 1 everywhere in the space-time 
if (am) « 12e. The last condition can always be held by an appropriate choice of 
the integrating constant (am). For N = 24 (c = 1/2) the curvature R was shown 
to be constant: R = e~( w+1 ^ . By an appropriate choice of the constant w one can 
control the condition R << N, so the semiclassical approximation is correct also in 
this case. 



1 



r 1 j_ + 
[- + x^ - x [ 



ol, if A^a 



(5.44) 



4em a 
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6 Discussion 



In resume, we have obtained that the preliminary hope that one-loop quantum 
corrections remove the classical black hole singularity is not realized for very large 
N (c >> 1). The space-time singularity is still present in the general solution for 
the action (5.2) in this regime. However, something interesting happens when N 
takes some finite (not very large) values. We have shown that for N = 48 (c = 1) 
the solution of the system (5.2) describes geodesically complete space-time regular 
everywhere. The corresponding scalar curvature (5.41) takes only finite values. In 
the other case, when N = 24 (c = 1/2) the semiclassical action (5.2) describes 
the (de Sitter) space-time of constant curvature. This space-time is also obviously 
free from singularities. Remember that N is the number of scalar fields or, more 
generally, N is the number of sorts of particles in a matter multiplet. 

We conclude with some remarks in the order of discussion. It seems to be 
reasonable to consider the requirement of space-time regularity as some kind of 
principle: " The space-time singularities must be absent in the complete quantum 
theory of gravity and matter " . Then, our semiclassical analysis can be interpreted 
as that this "regularity principle" is not valid in general. But it leads to some 
restrictions on the particle contents of the theory. In the case under consideration, 
it constraints the number of matter fields N. There are some well known principles in 
modern physics which bound the particles spectrum: for example, the requirement 
of anomalies cancellation. Therefore, it would not be very surprising if the black hole 
physics gives us one more. In this paper, we have considered the two-dimensional 
case. However, the same situation can take place in four dimensions [15]. 

Of course, our study is just semiclassical and cannot be considered as a strict 
proof. The analysis in the framework of the complete quantum theory is necessary. 
However, the above consideration seems to be an argument in favor of the hypothesis 
on the relation between absence of the space-time singularities and particle spectrum 
of the theory. 
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8 Figure Captions 



Fig.l: The shape of the function (f>(R) in the vicinity of point R = i?o where 
<fi'(R) = f"(R) = 0. There are two values R a , Rb of R which correspond to the same 
value 4>. So the inverse function R(<p) has two branches. 

Fig.2: The space-time near the time-like line R = Rq ((f) = 0) where (j>'(R) = 
f"(R) = 0. It consists of two sheets glued along the line R = Rq. 

Fig. 3: The space-time obtained by the falling of <5-like impulse of matter at 
x + = Xq on the black hole. For x + > x^ the singularity is slightly deformed and 
asymptotically reaches the new horizon at x~ = x^. 

Fig. 4: The black hole space-time deformed by quantum backreaction for c > 1. 
The singularity now is time-like and points of horizon satisfy the condition x + x~ = 
0. The regions I and III are asymptotically flat. 

Fig.5: The space-time obtained by the falling of <5-like impulse of matter at 
x + = Xq on the black hole for c > 1. The singularity for x + > x^ is slightly shifted 
and asymptotically tends to new horizon at x + = x~\ and x" = x^ . 
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